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ABSTRACT 



Context. The evidence of a decrease with increasing luminosity in the fraction of absorbed and Compton-thin among X-ray- 
selected (2-10 keV) AGN is observationally well supported, while that of an increase in /n^, with redshift is fairly controversial. In 
Lamastra, Perola & Matt (2006), the gravitational effect of the SMBH on the molecular interstellar gas, in the central region of the 
host galaxy, was shown to predict an anti-correlation between fats and the black-hole mass Mbh- 

Aims. The most recent findings on the distribution of the Eddington ratio X=Li,ILe as a function of Mg^ and z are used to convert that 
relationship into one between fabs and both bolometric {Lt) and X-ray (Lx) luminosities at various values of z- 
Methods. The findings for A(Mbh,z) are properly treated to ensure completeness in the prediction of /a;„ above a certain luminosity, 
at values of z=0.1, 0.35, 0.7, and >1. To verify the consequence of these findings alone, we first adopted a distribution of gas surface 
density E, observed in a sample of local spiral galaxies, irrespective of the galaxy morphological type and z- 

Results. Assuming the Eddington limit, A=l, in the A(Mbu,z) distribution as a "natural" cut-off, the predictions are consistent with 
the existence of an anti-correlation between fahs and Lx, but they fail to reproduce an increase in /„(,, with z- Because the early type 
galaxies are on average much poorer in molecular gas than late type ones, a quantitative agreement with the local value of f^bs requires 
the existence of a correlation between S and the central activity. An increase in typical values of S with z, correlated with the activity, 
might explain an increase in f^bs with z- However, f^bs could hardly exceed about 0.3 at the highest luminosities. 

Key words. Galaxies: active - X-rays: galaxies - ISM: clouds 



1. Introduction 

The selection of active galactic nuclei in the hard X-ray band (up 
to about 10 keV with the XMM-Newton and Chandra satellites), 
together with their redshift z measurement after optical identifi- 
cation, immediately yields the spectroscopic estimate of the pho- 
toelectric ally absorbing column density, Nh, intrinsic to each of 
^ them. A schematic split between objects with Nh either smaller 
• ^ or larger than 10^^ Hcm^^ corresponds fairly well (although not 
exactly one-to-one) with the optical classification of AGN as ei- 
H ther type 1 (broad-line) or type 2 (naiTow-line) (see, e. g., Perola 
. . . et al. 2004). Among the objects with Nh larger than 10^^ Hcm'^, 
a further subdivision is necessary, those with Nh smaller than 
cr^' = 1.5 X Kfi'^arT^, called Compton-thin, and those with Nh 
> cTj^ , called Compton-thick. The latter category is hardly repre- 
sented in samples extracted from images with maximum X-ray 
energy about 10 keV; therefore, the investigations on the frac- 
tion of absorbed AGN do in fact apply only to the Compton-thin 
sources: /„;,, = n{\Q^^ < Nh < lQ^'^)/n{total). 

The empirical results for this fraction, reported in the liter- 
ature, have not yet reached full agreement. To first mention the 
two papers (Ueda et al. 2003, La Franca et al. 2005) they ad- 
dressed the issue of fats in detail simultaneously with that of the 
AGN luminosity function in Lx{2 - 10 keV), Lx hereafter, and its 
cosmological evolution, they agree on the existence of an anti- 
correlation between fabs and Lx, but not on its behaviour with 
z. While its normalization and shape stay constant,, according to 
Ueda et al. (2003), according to La Franca et al. (2005, LF2005), 
fabs{Lx) increases with z and becomes shallower, so that the in- 
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crease is more marked at high than at low luminosities. Akylas 
et al. (2006) confirm the existence of the anti-correlation with 
Lx, but argue that the apparent increase in fabs with z in their 
data sample might be due to a systematic overestimate of Nh at 
high redshifts. Evidence of the anti-correlation with Lx has also 
been reported by Hasinger (2004) and by Treister & Urry (2005, 
see also Bargeret al. 2005). To avoid the uncertainties associated 
with the estimate of Nh, Treister & Urry (2006) use the optical 
spectroscopic classification of the counterparts in a combination 
of seven X-ray-selected samples and claim evidence (after cor- 
recting for selection bias) of an increase with z of the fraction of 
obscured AGN. They do not, however,explicitly address the is- 
sue of whether the shape of the anti-correlation might change as 
well. A correlation between the fraction of broad-line AGN and 
their luminosity is reported by Simpson (2005), after analysing 
an optically selected sample (z <0.3) and using the [OIII] narrow 
emission line to indicate the bolometric luminosity. 

Although Dwelly & Page (2006), after a detailed study of 
one X-ray deep and naiTow beam survey, find no evidence 
of both the anti-correlation and of a redshift dependence in 
fabsiLx), we consider it reasonable at present to adopt the view 
that the existence of the anti-correlation between fabs and Lx 
is fairly well-supported. For X-ray-selected samples in particu- 
lar, we regard the results as more trustworthy that are based on 
the combination of deep-narrow and shallow-wide surveys (see 
LF2005 for a discussion on this point). On the other hand, we re- 
gard the evidence of an increase in fabsiLx) as still controversial 
with increasing z, even more so that of a change in its "slope". 
In a previous paper (Lamastra, Perola & Matt 2006, Papl), we 
proposed a simple model to explain the anti-coiTelation, which 
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calls into play interstellar gas for the Compton-thin sources, 
while the molecular torus, as invoked in the standard "unifica- 
tion" model (Antonucci 1993), remains the best explanation for 
the Compton-thick sources. Our proposal is no more than an ex- 
tension of the standard model, whose best empirical support is 
the finding by MaioUno & Rieke (1995) of a correlation between 
absorption and disc inclination in (optically classified as) inter- 
mediate type Seyfert galaxies, which is absent in extreme type 2 
Seyferts. A finding whose simplest explanation is that the torus 
and the disc are, in general, not necessarily coplanar. The model 
assumes that the obscuring matter is distributed in a rotation- 
ally supported disk, extending out to a few hundred of parsec, as 
observed in nearby spiral galaxies mainly in the form of molec- 
ular hydrogen. The gravitational pull of the central supermas- 
sive black hole (SMBH), combined with that of the bulge (see 
Ferrarese & Ford 2005 for a review of the connections between 
these two galactic components) on the gaseous disc, as a func- 
tion of the distance R from the centre, reduces the gaseous disc 
covering factor C as the SMBH mass, Mbh, increases. This ef- 
fects leads to prediction of an anti-correlation between fabs and 
Mbh, which is likely to be reflected into the one observed be- 
tween fabs and Lx- In Papl it was shown that the latter could be 
reproduced fairly well in the "local" universe by adopting a "typ- 
ical" value of 0.1, irrespective of the mass, of the ratio between 
the bolometric and the Eddington luminosity, A = Lb/Ls. This 
value of A was consistent with the estimates available then (see 
references in Papl). 

Recently, by means of a wealth of spectroscopic data from 
the Sloan Digital Sky Survey (SDSS), the empirical knowledge 
of the behaviour of /I as a function of Mbh has had a substantial 
turn, thanks in particular to the work by Netzer & Trakhtenbrot 
(2007, NT2007) on broad-line AGN. They explored it for masses 
from 10^ to lO' Mq in the z interval 0.3-0.7, and found a strong 
dependence on both M/;// and z- It should be noted that several 
earlier papers seemed to agree on a practically constant typical 
value of A with increasing z, irrespective of luminosity and mass, 
after selection effects had been taken into account (a good ex- 
ample, over the range z=0.3-4 in Kollmeier et al. 2006, see also 
Lamastra, Matt & Perola 2006). In the framework of a model 
where, by definition, the intrinsic properties of the central en- 
gine are the same in broad and narrow-Une objects, it is nat- 
ural to assume that the behaviour found by NT2007 holds for 
the narrow-line AGN as well. As a matter of fact, in the "local" 
(around z=0. 1 ) universe, a qualitatively similar result for A(Mbh) 
was found earUer on by Heckman et al. (2004, H2004) in a large 
sample of SDSS selected narrow-line AGN. Despite some criti- 
cism of the reliability of the methods used by H2004 expressed 
by Netzer et al. (2006), we use their findings as well. 

The new light on the Mbh-A-z relationships in the range of z 
up to 0.7 prompted us to revisit the model of Papl, with the fur- 
ther aim of verifying whether the z-dependence of fabs claimed 
by LF2005 might be reproduced or otherwise. 

This paper is organised as follows. Section 2 is devoted to the 
new results mentioned above and to their quantitative represen- 
tation as is used later. Sect. 3 to the mass function of the SMBH, 
Sect. 4 to the molecular gas surface density in the inner region 
of galaxies, both late and early-types. Sections 5 and 6 present 
the model "predictions" of fabs as a function of Lboi and of Lx, 
and their comparison with the LF2005 results. A discussion and 
a conclusion follow in Sects. 7 and 8. 



2. New light on the Mbh-A-z relationships 

NT2007 uses a sample of almost ten thousand broad-line AGN, 
with m, < 19.1, from the SDSS to investigate the distribution of 
A in different range values of Mbh, from 10^ to lO'M©, as a func- 
tion of the redshift z- They estimate Mbh using the L5100 - Rblr 
(the radius of the broad-line region) relationship (as given by 
Kaspi et al. 2005) in combination with the FWHM(//j8). The 
value of A is obtained after applying a (constant) bolometric cor- 
rection to Lsioo- The results, as summarised in their Table 1, ap- 
ply to intervals of z up to Zmax=Q-1^, except for the mass range 
around Iff Mq for which Zmax=0-^, and down to Zm,„=0.05 for the 
mass ranges around 10^ and IO^^Mq, with Zmin increasing up to 
0.25 for the mass range around IO'Mq. Thus, in the z interval 
0.25-0.75, file results apply to all masses from lO''-^ to IO'Mq, 
while the overlap is confined to the narrow z interval 0.25-0.3 for 
the mass range around 10^ Mq. These results concern a feature, 
'^peak, defined as the maximum in the observed ^-distributions, 
and its dependence on Mbh and on z. NT2007 apply reliable sta- 
tistical tests to demonstrate that the dependence found is most 
Hkely real, rather than a consequence of the flux Umit of the 
sample. 

To avoid misunderstandings on the extent to which we rely 
on NT2007 in the following , it is appropriate to comment on 
the "incompleteness" of their sample. As extensively discussed 
in particular by Richards et al. (2006), whose aim was to esti- 
mate the luminosity function of QSO as a function of the red- 
shift, samples of such objects from the SDSS are affected by 
incompleteness due to the target selection algorithm, based on 
the observed photometric colors, and to the limiting magnitude. 
From the statements in NT2007, we infer that the value of A 
(see the marks in their Fig. 3, call it Amark), below which they 
regard their sample incomplete, relates to the limiting magni- 
tude OT, < 19.1. Figure 6 in Richards et al. (2006) shows clearly 
that the incompleteness degree increases very steeply for mag- 
nitudes above this limit. In the same figure it is apparent that 
the incompleteness due to the color selection is of the order of 
20% (averaged over 15< m,- <19.1) and stays constant up to z 
about 2. Going through the way Richards et al. (2006) then used 
these results to obtain the proper normalization of the luminosity 
function, we are led to conclude that the color-induced selection 
effects do impact on the "absolute" values, not on the "relative" 
ones. In other words, the "shape" of the distributions of A above 
Amark should uot be affected by these effects. Since for our pur- 
pose, as will be apparent later, it is the "shape" of these distribu- 
tions that is essential, we now describe the way we have adopted 
to quantify this shape. 

We introduce the distribution function <I>{Mbh,A,z)- Figure 
3 in NT2007 displays the fraction of AGN as a function of A 
for two ranges of black hole mass (lO^-^-lO^ ^ lO^-^-lO^-^ Mg), 
each one for three values of z. Except for the range MBH=iO^'^- 
IqT-^ Mq at z=0.5, A„u,rk is placed somewhat to the left of the 
peak in the "observed" distribution in the other five cases. We 
thence looked for an analytical function that could represent the 
five "observed" distributions above Amark reasonably well. We 
eventually adopted a Lorentzian profile: 

<^{Mbh, a, z) = ^ Jjrj^ A > ^mark, (1) 

with a "natural" cut-off at /i=l. 

After assigning statistical error bars corresponding to a total 
number of 150 objects (a relatively arbitrary choice between the 
total numbers from 50 to 230 declared in NT2007 without fur- 
ther specifications)to each one of the five histograms in NT2007, 
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Fig. 3, we operated ax^ fit to obtain the five values of Apeak and 
of w. After noting that all five values of w turned out to be equal 
to Apeak multiplied by approximately the same constant, we pro- 
ceeded with a simultaneous fit after requiring that this constant 
should be exactly the same. The result obtained is w = 1.5 Apeak, 
with a reduced ;^'^=0. 68. We then adopted this result to w to hold 
for any mass value from 10^ to 10^ Mq, independent of z- It is 
to be noted, at last, that our values of Apeak do not differ sig- 
nificantly from those, apparently obtained by NT2007 directly 
from the position of the peak in their fractional distributions, 
displayed in their Fig. 4. (We later return on the normalization 
A{Mbh,z).) 

For the z-dependence of ApeakiMsH), NT2007 find (see also 
their Fig. 4) that it can be represented (of course only in the 
intervals of z effectively covered) with a power law, namely: 

Apeak ^ Z'"^""'^ (2) 

where the exponent y increases with Mbh- We therefore adopt 
the following z-dependence: 

log ApeakiMBH, z) = log a(MBH) + tC^bh) log z. (3) 

For our application, it is convenient to assign an analytical 
form to the function y{MBH) and aiMBH)- For the former, it turns 
out that the five values of y in NT2007, Table 1, are well-fitted 
with 




log z 

Fig. 1. The redshift dependence of Apeak for five values of Mbh in 
Mq, derived from NT2007 as described in the text. The crosses 
are the values obtained by fitting the Lorentzian profile ([T]l, the 
hexagons are from Fig. 4 in NT2007 



yiMBH) = -3.01 + 0.581og(MB,f). 



(4) 



For aiMBH), which is not given in NT2007, we have to rely upon 
the only two values obtained with the above-mentioned fits. We 
choose to connect them with a straight line in the log-log plot of 
a and Mbh, and obtained 



a{MBH) = IQ'\MbhIMJ^-'\ 



(5) 



The power-law relationship (O is illustrated in Fig. [T] for five 
representative mass values, from 10^ to 10^ Mq. Except for 
Mbh=^Q^ Mq, the normalization comes from (|5]l. The extrap- 
olation of Eq. ^ to this mass places the power law well above 
the three Apeakiz) points in NT2007, Fig. 4, which are reported 
in our Fig. [T] We therefore decided to adopt a power law with 
a slope as given by Eq. (01 and a normalization in agreement 
with these points. The figure illustrates very well that the power 
law (O pertaining to a given mass always stays above the one 
pertaining to the next higher mass. 

To use the function there are three other elements to be 
taken into consideration. The first would be the extrapolation 
of below Apeak{MBH,z). This is extremely uncertain, so we 
avoided even trying. We exploit the fortunate circumstance that 
the information available on the shape of O above Apeak will al- 
low us (see Sects. 5 and 6) to predict the quantity fabs for L^o/ > 
IQ/^^ergs^^ (corresponding to ^ IQI^^-^ergs^'^ , Marconi et al. 
2004), namely the range of where the increase in fabs with z,is 
most pronounced according to LF2005. Incidentally, this is also 
the luminosity interval where the SDSS quasar luminosity func- 
tion was determined by Richards et al. (2006); see also Hopkins 
et al. (2007). 

The second is the normalization A(Mb//,z). Since we use O 
only for A > Apeak, it is convenient to introduce the integral quan- 
tity 



^(Mbh, a, z)dA = q(MBH, z) 



(6) 



This quantity, which depends on A(Mbh, z), represents the frac- 
tion of SMBH, as a function of their mass and cosmological 
epoch, which are active above Apeak- 

There are no elements in NT2007 to infer the value of q. For 
our purpose, however, what matters is whether this fraction is 
either a constant or a function of Mbh', in other words, whether 
it preserves the same shape as the active SMBH mass function 
[q = constant) or not (q depends on Mbh)- In Sect. 5 we illustrate 
the consequences of both types of assumption. 

The third is the mass function (MF) of the SMBH in galaxies, 
F(Mbh), in particular that of the active ones. The MF is needed 
to weight the fractional function O when the contributions to 
fabs by different masses are combined within the same bin of 
luminosity. The MF will be dealt with in Sect. 3. 

We lack a similar analysis for A{Mbh, z) carried out explic- 
itly on broad-line AGN in the low-redshift universe. Thence we 
resorted to the work by H2004 on a large (23,000) sample of 
narrow-line AGN drawn from the SDSS. To the extent that this 
other type of AGN can be assumed to share the same intrinsic 
properties of the broad-line type, we can confidently use their 
results to include the local (out to z=0.1) universe in our work. 

Unlike NT2007, H2004 estimate the mass using the Mbh - 
cr* (stellar velocity dispersion in the galaxy bulge) relationship 
(as given by Tremaine et al. 2002). The need to measure cr* prac- 
tically excluded from the sample the brightest Seyfert galaxies 
and the QSO, a cautionary remark emphasised by H2004. They 
derive Lboi from Lgni multiplied by a constant factor (3,500, es- 
timated variance of 0.38 dex). This procedure is not completely 
reliable, because there seems to be a correlation between the cor- 
rection factor and the luminosity (Netzer et al. 2006), but we re- 
gard it as acceptable in the first approximation for our purpose. 

The quantitative results that we use are those summarised 
in H2004, Fig. 3 (left panel), where the integral distribution in 
A of the fraction of SMBH, which are active, is given for six 
values of Mbh, from 3x10'' to lO'' Mq. It can be immediately 
seen that, at least qualitatively, the trend of A with Mbh is the 
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same as found by NT2007. This form of presentation, namely 
the fraction of active objects with respect to the "local" F{Mbh), 
is particularly convenient because it can be directly used, as we 
do in Sect. 6, without resorting to assumptions on the quantity 
q, Eq. (|6]l. It is also interesting to check whether these distribu- 
tions can be represented by the function ([T]l- Since the curves in 
H2004, Fig. 3, are devoid of error bars, we could only adjust, 
as best as we could, the integral of this function to them. The 
outcomes are in general quite satisfactory, and one example is 
shown in Fig. |2](where the curve from H2004 is represented by 
a selection of points). Typically the high A tails stay somewhat 
above the H2004 curves. One might require that the excess in the 
Lorentzian tail does to some extent take care of the lack of the 
brightest objects in the original sample, but we need not insist 
on this issue, because the difference is of little consequence for 
our final results. 




-2.5 -2 -1.5 -1 -0.5 
log X 



Fig. 2. Integral distribution in A of the active fraction of SMBH, 
Mbh=10^ Mq, out to z=0. 1 . The crosses are derived from H2004, 
Fig. 3, the line is our Lorentzian representation. 



The parameters extracted from H2004 that we use in Sect. 6 
are summarised in Table [T] The Lboifi in the second column (with 
the corresponding /Iq in the third) represents a completeness 
limit valid for all masses between 10^ and 10^ M©. The fourth 
and fifth columns contain, repectively, A,„ax, the maximum value 
that can be read in H2004, Fig. 3, and qn, the fraction of active 
SMBH above Aq from the same figure. From the sixth column 
onward the quantities given were obtained using the Lorentzian 
profile, as explained above. The quantity Apeak is very low and 
badly determined, but the quantity that really matters is w; the 
fraction qni in the last column differs only marginally from qn- 

It should be noted that the fraction of active SMBH in H2004 
is given with respect to the local F(Mbh), such as the one given 
by Marconi et al. (2004) and described in the next section. 

3. The mass function of SMBH 

The number density of SMBH as a function of mass, F(Mbh), 
has been estimated in the local universe by several authors (see 



Marconi et al. 2004, Shankar et al. 2004 and references therein). 
The one adopted here is from Marconi et al. (2004) and is shown 
in Fig. [3] where the MF which includes all morphological types 
(their Fig. 2b) is given along with the MF in late, Fl, and early, 
Fe, type galaxies separately (A. Marconi, priv. comm.). This Fl 
is preferred to the one estimated by Shankar et al. (2004) for the 
attention paid to the change in the bulge to the total luminosity 
ratio along the Hubble sequence of late type galaxies. The reason 
we keep the late separate from the early type galaxies will be 
explained in Sect. 7. 

For our purpose, what does matter is the shape of the MF. 
Insofar as the growth of the SMBH is a consequence of accretion 
and is therefore accompanied by the AGN type activity, it is most 
likely that the shape of the MF for the "active" SMBH changes 
with redshift and is different from that of the whole population 
as observed in the local universe. In principle, the redshift depen- 
dence of the "active" SMBH mass function, F*(Mbh,z), could 
be inferred from the AGN luminosity function as a function of z, 
provided that a simple relationship holds between Mbh and the 
luminosity (for instance, A = constant). In practice, as discussed 
in Marconi et al. (2004) and as has become particularly evident 
after NT2007, the conversion is neither straightforward nor uni- 
vocal. Hence F*{Mbh) can only be reliably estimated through 
a "direct" evaluation of Mbh in a properly selected sample of 
AGN, with special care devoted to correcting for incomplete- 
ness selection effects. On the basis of SDSS samples, H2004 
(z < 0. 1) for naiTow emission line AGN and Greene & Ho (2007, 
z < 0.3, with masses derived as in NT2007) for broad emis- 
sion line AGN, both find that in the local universe F*(Mbh) is 
steeper, above 10^ Mq, than F(Mbh)'- a difference attributed to 
the so-called downsizing 

In Sects. 5 and 6 we, in a first instance, adopt the "shape" 
of the local F(Mbh)\ however, the impact on the results of other 
options for the downsizing is also discussed. 




6 7 8 9 10 
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Fig. 3. The local SMBH mass function from Marconi et al. 
(2004). Dotted line: late type galaxies; broken line: early type 
galaxies. 
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4. The molecular gas surface density £ 

The basic feature of the model in Papl is the gravitational force 
exerted by the BH (and by the associated stellar bulge) on a ro- 
tationally supported disc made of molecular clouds. This force 
shapes the gas distribution profile, as a function of the radial dis- 
tance R from the BH, in such a way as to determine an anticorre- 
lation between Mbh and the "covering factor" determined by the 
disc. Quantitatively this factor was calculated in Papl assuming 
that the distance out to which the disc extends is at least equal 
to 2Ri„fi (Rinfi is the distance from the BH at which the molec- 
ular gas distribution profile has a point of inflection, see Eq. (9) 
in Papl), that is, about 25 pc when MBH=iO^ and extending to 
about 450 pc when Mbh=10^ Mq. 

Figure|4]shows this "covering factor", C(Z, Mbh), as a func- 
tion of the surface density S, from to 2000 Mq pc^^, for four 
values of Mbh, from 10^ to 10"^ Mq. After a steep growth, the 
curves enter a regime where the increase is very slow. The rather 
sharp change in slope takes place around a S value from about 
200 to about 400 Mq pc^^ when moving down from Mb//=10' 
to Mbh=^0^ Mq. Thus, in predicting fabs, attention must be paid 
to the distribution of E adopted, particularly on the fraction of 
objects below those values. 




2000 



Fig. 4. The covering factor C as a function of S, according to 
the model, for different values of the SMBH mass. From top to 
bottom: 10^ 1Q\ 10^ and lO'' Mq. 



In Papl (to which we refer the reader for more details), we 
used the only results obtained with sufficient angular resolution 
in a molecular line. These are represented by a sample of 44 



nearby and bright spiral galaxies mapped at the 3 mm CO J=l- 
line, the BIMA SONG survey (Heifer et al. 2003). At their 
average distance of 12 Mpc, the angular resolution of 6 arcsec 
corresponds to a radial distance from the centre of about 210 pc. 



1000 



100 r 




Fig. 5. The molecular gas central, face-on surface density S of 
the spiral galaxies in the BIMA SONG survey, versus the mass 
estimate of the SMBH hosted by these galaxies. Hexagons: Mbh 
from the Mbh - Lbuige relationship; open: from Dong & De 
Robertis (2006); filled: our work. Triangles: from the Mbh - o"* 
relationship, Batcheldor et al. (2005), Merloni et al. (2003), Woo 
& Urry (2002), Pellegrini (2005). Crosses and pentagon: from 
maser and from stellar kinematics. Ho (2002). Square: from re- 
verberation mapping, Kaspi et al. (2000). 



In Fig.|5]the central face-on S of these objects is plotted as a 
function of the mass of the SMBH (estimated as explained in the 
figure caption). For the majority of these galaxies Mbh is com- 
prised between 10^^ and 10** Mq. In practice this is the range 
of masses for which the information described in Sect. 2 will be 
used to predict fabs in Sects. 5 and 6. Given the small number of 
objects in this sample, we use the observed distribution of E in 
the 44 objects to evaluate the uncertainty on fabs- 

In Papl we stressed that the effect under study does not re- 
quire any type of correlation between the central activity and 
the presence of the absorbing gas, the latter being "in place" as 
a property of the galactic interstellar medium. This statement 
can be regarded as basically correct for the late type galaxies 
(not only in the local universe, as shown by the survey just men- 
tioned, but presumably also at cosmological distances), and now 
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we need to consider to what extent it might be applied to the 
early type galaxies. 

Elliptical (E) and lenticular (SO), unlike the spiral (S) galax- 
ies, are rare in the field and frequent (the bright E's in particular) 
in galaxy groups and clusters. Therefore, useful interferometric 
observations are only available for a dozen of them, optically lu- 
minous or very luminous (thus hosting SMBH well above 10^ 
Mq), out to a maximum distance of about 85 Mpc and with a 
strong CO line signal. It is convenient to concentrate first on 
these results, and we limit ourselves to those obtained by Young 
(2002, 2005) with the BIMA interferometer on seven "true" E 
(that is galaxies with an r^^'^ photometric profile). Molecular gas 
is found distributed in rotating, mostly symmetric discs, with 
radii ranging from 1 to 6 kpc. For each of them we report the 
peak surface brightness, translated into a surface density of H2 
gas, along with the galaxy surface area corresponding to the 
beam. UGC1503 (field): 60 pc^^ (2.4x2.1 kpc^); NGC807 
(field): 40 Mq pc^^ (2.4x2.0 kpc^); NGC3656 (merger rem- 
nant): 720 Mo pc-2 (1.7x1.3 kpc^); NGC4476 (Virgo cluster): 
112 Mo pc-2 (0.73x0.49 kpc"); NGC5666 (field): 115 Mq pc^^ 
(1.7x1.3 kpc2) (from Young 2002). NGC83 (group): 135 Mo 
pc-2 (2.7x2.5 kpc^); NGC2320 (A569): 353 Mo pc^^ (2.4x2.2 
kpc^) (from Young 2005). On account of the values of the 
area/beam, it is evident that these selected cases correspond to 
the values of S, which imply a substantial covering factor C and 
which occur in a large fraction of the S galaxies in Fig. |5] But 
how frequent the occurrence is of such a situation in E and SO 
galaxies? The answer to this question does not appear to be well 
settled yet. Furthermore, most data are from single-dish obser- 
vations, which can be properly converted only into a total H2 
mass. 

On this issue it seems to us appropriate to first quote the com- 
pilation by Bettoni, Galletta & Garcia-Burillo (2003) of a wealth 
of data obtained prior to the year 2003. Most useful for our pur- 
pose are this mean sample values, subdivided per morphological 
class according to the mean numerical Hubble stage index t as 
tabulated in the Third Reference Catalogue of Bright Galaxies 
(RC3 de Vaucouleurs et al. 1991), of the quantity M(//2)/Z)^g, 
where D25 is the linear size corresponding to the 25 mag foresee^ 
isophote. After assembling their classes into three groups, Gl 
(f from -5 to -3, E), G2 (f from -2 to 0, SO), and G3 {t from 
1 to 6, Sa to Sc), from Table 5 in Bettoni, Galletta & Garcia- 
Burillo (2003), we obtain the following proportions for the quan- 
tity M(//2)/Z)2g, G1:G2:G3 = 0.06:0.29:1. Thus, with the cau- 
tionary remark that the size of the H2 gas structure in the three 
groups may not scale similarly with D25, the immediate conclu- 
sion one can draw is that in the ellipticals the surface density of 
molecular gas is, on average, more than one order of magnitude 
less than in spirals, while it is a factor of about three less in the 
lenticulars. This conclusion we believe is stronger than the one 
that can be drawn from the detection rate, which is very sen- 
sitive to the sample selection, as discussed in Combes, Young 
& Bureau (2007); see also Welch & Sage (2003), Sage, Welch 
& Young. (2007) and references in these papers. An additional 
important result emerging from these papers is that, among the 
early type galaxies, the ratio M(//2)/i' decreases with increasing 
luminosity L, despite the increase in depth of the potential well. 
Taken together with the results quoted from Bettoni, Galletta & 
Garcia-Burillo (2003), this finding leads us to conclude that the 
more luminous (and therefore hosting the more massive SMBH) 
an early type galaxy, the less the chance of finding "in place" H2 
surface densities similar to the ones reported above from Young 
(2002, 2005). 



For high z objects, we would like to mention the most obvi- 
ous argument used to account for the Mbh - c* and Mbh - L 
observational relationships (see Ferrarese & Ford 2005 for a re- 
view), where cr* and L are the stellar velocity dispersion and 
the luminosity, respectively, of galactic "bulges" or "spheroids". 
The argument goes that the formation of the spheroids is closely 
connected with the growth of the SMBH sitting at their cen- 
ters. If the two phenomena had progressed hand-in-hand from 
the earliest to the present times, it would be legitimate to pos- 
tulate that a phase of intense activity (AGN phase) should in 
general be accompanied by star formation on the spatial scale 
of the bulge. The latter requires the accumulation on this scale 
of gas (whatever its origin) in a dense, most likely molecular 
state, wherein the star formation takes place. From an obser- 
vational point of view the validity of such a postulate rests on 
results obtained with studies of vast samples of galaxies, based 
on SDSS spectrography and photometry. Here we refer in par- 
ticular to Kauffmann et al. (2003a), Kauffmann et al. (2003b), 
H2004. Quoting from the review by Heckman & Kauffmann 
(2006): "Strong emission-line AGN inhabit those unusual galax- 
ies that are both relatively massive and dense, yet have a sig- 
nificant young stellar population". The denomination "unusual" 
emphasises the contrast with the usual make-up of galaxies with 
similar mass and structure, which, within the limits in sharpness 
of the observational distinctions employed, contain only a popu- 
lation of old stars. 

A cautious theoretical elaboration on these results must in- 
evitably contemplate the possibility that the SMBH activity is 
a consequence of the transformation into stars of the dense gas 
and the dynamical implications of supernovae explosions in par- 
ticular. In this instance, it could well be that, during most of the 
active phase of the SMBH, the dense gas is no longer present 
in sufficient quantity on the spatial scale required by our model. 
Furthermore, it is doubtful that the structure of the molecular gas 
were dominated, in such a phase, even if present, by a rotation- 
ally supported configuration out to several 100 pc. Only direct 
observations in molecular lines, with the adequate angular reso- 
lution, will answer these questions. 

To reach simple and sharp conclusions, we therefore con- 
sider two instances. In Sects. 5 and 6, in order to emphasise only 
the consequences of what was described in Sect. 2, we shall posit 
that, for all AGN no matter the morphological type of the host 
galaxies, the model applies equally well and a single represen- 
tative distribution of 2 will be used. In Sect. 7, we posit that the 
above applies only to the late type galaxies, while in the active 
early type galaxies the "covering factor" is basically unknown. 
The latter instance requires that Fi{Mbh) be kept separate from 

FEiMsH)- 

5. Predicting as a function of iuminosity for 

Mfiff > 10' Mq, at z=0.35 

The fraction of absorbed AGN, fats, in a given luminosity bin 
(Lii < L,!, < Lfe2) can be computed using the equation 

T.iF{MBH,d^MBH,i ^^y{MBH,i,A.,z)dA 

^here = t^^^) (^1 > V-*) and ^2 = T^zmin^ ^ 
/(Ej) is the fractional distribution in the quantity E, introduced 
to take properly into account the strong dependence of C(2) on 
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low values of the argument. In our application we used mass bins 
0.3 dex wide, and the results are presented in bins of L/, that are 
0.5 dex wide. A normalization of the function O, as explained in 
Sect. 2, will depend on the assumption we make on the quantity 
qiMsH) given by Eq.®. 

To predict fahsiLt) with Eq. (|7]l, we only consider the ac- 
tive SMBH with A > Ap^akiMBn)- When the relationshps (3) and 
(5) are combined, it turns out that, at one particular value of the 
redshift, z=0.35, ApeakiMBH) °^ ^bh ^^^y closely. Therefore, all 
masses in the range 10^'^- 10^ M©, to which at z=0.35 both (3) 
and (5) apply, contribute above the same value of Lh, equal to 
\(f^ergs~^ . In other words, the active BH in this range of mass 
with A < Apeak contributes only below lO^^ergi"'. Notably, this 
is the Ledd of Mbh = 10^ Mq. Thus, to obtain a "complete" es- 
timate of fahs in the interval of logL/, = 45-45.5, only the con- 
tribution of the masses in the range 10^-10^'^ M© needs to be 
added, which we did by resorting to a small extrapolation of the 
curve labelled 10^ Mq in Fig. 1 beyond its formal limit at z=0.3. 
This convenient situation will be exploited to make some clear 
point, and to predict fabs at a given z to be used in Sect. 6 for 
comparison with the predictions at higher and lower values of 
the redshift. 

As anticipated in Sects. 3 and 4, we first adopt the "shape" 
of the local F(Mbh), Fig- S as representative of the ensamble 
of active SMBH in both late and early type galaxies, and assume 
that the "model" applies equally well to the two types of hosts 
with the same distribution of S. Concerning the quantity q (Eq. 
Q), we start with the assumption that it is equal to a constant, 
namely independent of Mbh- It is perhaps worth emphasising 
that the value of this constant is irrelevant for obtaining the frac- 
tion fahs, while its precise knowledge would be necessary if we 
wanted to predict the number, for instance per unit volume, of 
the two types of AGN. The result is given as a histogram in Fig. 
|6l The first aspect deserving a comment is the sensitivity of the 
result to the choice concerning E. The values of fabALh) in the 
individual bins of the histogram have been computed using the 
"observed", local, distribution of the 44 values of E referred to in 
Sect. 4, and the bars represent the 68% confidence intervals cal- 
culated on the basis of the limited number of objects in that dis- 
tribution. For the sake of neatness, in the following similar his- 
tograms, the bars will no longer be reproduced. For comparison, 
in the same figure fabs is also given for two choices of a "single" 
representative value of Z, namely the median (150 Mq pc"^) and 
the mean (270 Mq pc"^) of the same distribution. In both cases, 
especially in the second, the prediction is significantly higher in 
all bins. This follows from the fact that (see Fig. HI, when a sub- 
stantial fraction of objects span the interval from to about 300 
Mq pc"^, the covering factor is extremely sensitive to the way 
the values of E are distributed. In this respect, the prediction of 
fabs obtained adopting E=270 Mq pc"^ (as it can be judged from 
Fig. |4]l is rather close to the ceiling attainable by pushing the 
distribution toward much higher values of E than adopted here: 
a point on which we return in Sect. 7. 

It is important to note that, without the "physical" limit 
of unity for A, the result would be different. For instance, if 
no cut-off at A=\ in Eq. ([T]i is imposed, the computation of 
fabs results is approximately a constant value, iiTespective of 
luminosity. This occurs because a single mass range, lO^-lO^ -' 
Mq, turns out to dominate at all luminosities, and the value 
of fabs to be close to the one pertinent to these mass values. 
Moreover, in this instance, we have no information to adopt 
from NT2007 in order to estimate the contribution of the 
masses below 10^ Mq, a very serious cause of uncertainties. 
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Fig. 6. The predicted fabs at z=0.35, obtained using the local 
SMBH MF, ^=const, with three different choices for E. Full line 
histogram: the observed distribution of E described in the text, 
the bars represent the systematic uncertainty due to the limited 
number of objects in this distribution. Dotted and broken line 
histograms: E=150 and 270 Mq pc"^ respectively. 



due to the relatively large number of SMBH in the mass range 
10''-10^ Mq. Because the relationship between Apeak and Mbh 
remains qualitatively the same at the redshift values dealt 
with in the following section, the same type of result, with the 
same type of uncertainties, is also obtained at the other redshifts. 

To illustrate the consequence of taking into account the 
downsizing effect mentioned in Sect. 3, we only need to change 
the leading assumption on q, namely adopt a proper dependence 
of q on Mbh- The shape of the broad-line AGN mass function 
in the local universe, as obtained by Greene & Ho (2007) above 
10^ Mq, suggests a dependence of the type q oc (M^/f/lO^)""'^. 
To show that the knowledge of the exact difference between F* 
and F is not very important for this range of mass, we assumed 
an even steeper correction as an exercise, namely q oc M^]^. The 
result is shown in Fig.|2] The anticorrelation becomes somewhat 
steeper, but the difference overall is not particularly marked. 
This is so because, even though the absolute number of objects as 
a function of Lb has changed radically, the fraction of absorbed 
sources depends solely on the mix of values of Mbh contributing 
to each bin of luminosity. In particular, for instance, in the high- 
est lumonosity bin the change is almost imperceptible because, 
despite the greatest change (decrease) in the absolute number of 
objects, the interval of Mbh values contributing to it have re- 
mained almost unchanged. 

To proceed to a comparison with the results obtained in the 
hard X-rays, we need to turn from Lb to Lj- The bolometric cor- 
rection used comes from Marconi et al. (2004) and is not linear 
We therefore recalculated Eq. (|7]) after converting Lb into L^. 
In the right panel of Fig. |7] the expectation for the two options 
on q is compared with the results obtained by LF2005 around z 
= 0.35. The two data points, with Poissonian error bars, repre- 
sent the objects in their samples within the interval z = 0.2-0.5, 
while the solid line represents the outcome at z=0.35 of their 
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Fig. 7. Left panel: full line histogram, same as in Fig.|6l dotted 
line histogram, q oc Mg\j. Right panel: histograms as in left panel 
as a function of Lx; solid line from the global fit by LF2005; the 
data points, with Poissonian uncertainties, include the objects 
with z=0.2-0.5 in the sample used by LF2005. 



"global" fit of the intrinsic distribution after correcting for se- 
lection effects. It must be noted that the "global" fit given by 
LF2005 (dotted lines in their figure 11) also includes Compton- 
thick AGN with 10^"* < Nh < 10^\ so we recalculated their 
best-fit relation to take only the contribution of the Compton- 
thin AGN into account. 

The difference between the data points and the line illustrates 
the relevance of the various selection effects that LF2005 take 
into account As it can be seen, our model predicts a decrease 
in the absorbed AGN fraction with increasing X-ray luminosity, 
but the slope of the predicted anticorrelation is not as steep as in 
LF2005. 



6. Redshift dependence in fabs 

We next move up to z=Q-T, the highest redshift where the 
NT2007 findings do apply. The results on fahs (Eq. U}, keeping 
q = constant, are shown in the left panel of Fig. [8] as a function 
of Lh, for Li, > 10"*^ '*erg5 where the completeness is ensured, 
and in the right panel of Fig. |8] as a function of L^. Notably the 
values of fabs are almost identical to those calculated at z=0.35 
for the same assumption on q. Despite the significant differences 
in ApeakiMBH) at the two redshifts, the parameters into play con- 
jure to give essentially the same quantitative outcome: not sur- 
prisingly in the highest bin of L^, which remains dominated by 
the highest masses. Even if, at this particular redshift, the global 
best fit at z=0.7 from LF2005 (solid line) appears to be in rea- 
sonably good agreement with our prediction, the fundamental 
point is that the z-dependence of the A{Mbh) distribution does 
not in itself introduce any increasing trend in fats, at least be- 
tween z=0.35 and z=0.7. This conclusion is strengthened by the 
fact that a possible dependence of q on Mbh, reflecting changes 
in the shape of F*{Mbh) above 10^ Mq, is likely to be less steep 
at z=0.7 than assumed at z=0.35, for demonstrative purpose, in 
the previous section, and would therefore make no significant 
difference for the conclusion. 

We now move to the local Universe. We used the informa- 
tion extracted from H2004, as described in Sect. 2 and sum- 
marised in Table [1] Figure |9] left panel, shows the results 
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Fig. 8. Same as Fig.|2l but only with ^=const, at z=0.7. The data 
points include the objects with z=0.55-0.85 in the sample used 
by LF2005. 



on fabs(Lb) obtained in two ways, one directly from the curves 
in Fig. 3 of H2004, the other from their Lorentzian represen- 
tation. The latter allowed us to go beyond logLi=46, although 
one should be aware that the absolute number of very lumi- 
nous objects is in any case very small in the local universe. Thus 
for our purpose, we consider uniportant that the H2004 sample 
lacks, by construction, the most luminous objects. In the interval 
logL/,=44.8-46, the fabs expected from the H2004 curves or from 
their Lorentzian representation are practically indistinguishable. 
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Fig. 9. Same as Fig.|7]atz=0.1, obtained using the H2004 results, 
see text. The dotted line histogram refer to the H2004 plots, the 
full line one to their representation with the function given by 
([Til. The data points include the objects with z<0.2 in the sample 
used byLF2005. 



After correcting Lb to Lx, the expected fabs is diplayed in 
Fig. m right panel. Here again the data points (z<0.2) and the 
best fit result from LF2005 are also given. The comparison 
shows that the predicted slope of the anticorrelation is definitely 
flatter than obtained by LF2005 at z=O.L 

To summarise these results, the fabs predicted as a function 
of Lb and Lx, for the three values of z covered (with ^=const for 
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z=0.35 and 0.7) are collected in Fig.fTO] The simple conclusion 
that can be drawn immediately is that, despite the strong depen- 
dence of A on Mbh and on z, the normalization and slope of the 
anticorrelation between f„bs and Lx, in the range of luminosities 
where they appear to change most quickly with z according to 
LF2005, is expected to change very little, if at all. 
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Fig. 10. The predicted fabs at z=0. 1 and, with ^=const, at z=0.35, 
0.7. In the right panel the dotted, solid, and dashed lines are from 
the global fit by LF2005 at z=0. 1, 0.35, and z=0.7 respectively. 



We conclude this section by tentatively pushing our predic- 
tions beyond z=l. In this respect we note that, up to z=0.7 (see 
Sect. 2) the A- Mbh relationships for the different values of Mbh 
do converge. Thus we see no contradiction with results, in partic- 
ular those obtained by KoUmeier et al. (2006) and mentioned in 
Sect. 1, going well beyond z=l, namely that out to z about 4 the 
mean value of A seems to stabilize at a value, <A>=0.25 indepen- 
dent of both Mbh and z, and that its distribution is approximately 
lognormal. A recent investigation by Netzer et al. (2007) does es- 
sentially confirm those results for what concerns the mean value 
of A in the redshift bin 2.3-3.4 and the mass range 10*"^-10'"^ 
Mq, while the distribution is somewhat broader than found by 
Kollmeier et al. (2006). A simplified approach, which ignores 
the distibution and associates the same value of A, equal to 0.25, 
to any value of Mbh > 10^ Mq, leads, when adopting the same 
2 distribution used previously, to the prediction of fabs given in 
Fig. [TT] which, by construction, is independent of the SMBH 
mass function, and hence of z- The only, albeit small, difference 
relative to the prediction at z=0.35, shown for comparison in the 
same figure, is that the expected fabs becomes somewhat steeper, 
contrary to the LF2005 results. 

7. Discussion 

Despite the uncertainty on the actual form of F*(Mbh,z) com- 
pared to that of the local F(Mbh), we have demonstrated that the 
predicted behaviour of fats at high luminosities does not repro- 
duce the increase with z found by LF2005, as long as the dis- 
tribution of E is kept fixed and indifferent to the morphological 
type of the host galaxy. This result can be understood easily, be- 
cause the trend in A{Mbh,z), in particular the way this quantity 
for different masses converge upwards with increasing z, implies 
that the higher the redshift, the higher the relative contribution 
of the more massive SMBH (with the lower C) at high luminosi- 
ties. Another way to appreciate this point is to note that, being 
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Fig. 11. The curve represent f„bs at z greater than 1, predicted af- 
ter adopting a fixed value of A=Q.25 irrespective of the SMBH 
mass. The histogram of the prediction at z=0.35, ^=const, is 
shown for comparison. 



Le(10^Ms)=10'^^ ergs-^ (Lx = lO'*'^-^ ergi"'), beyond this lumi- 
nosity only SMBH with higher masses contribute whose C stays 
in the range about 0.3; therefore, an increase in fabs above this 
range cannot be predicted. 

We are therefore left to operate on the other parameter of 
the model, namely E. On this quantity the empirical information 
we would need is practically missing, and only ad hoc assump- 
tions can be made. It should be noted immediately that, as shown 
in Fig. |6l the increase in a "typical" value of S shifts fabs up- 
wards, without markedly altering its slope, in fact just rendering 
it somewhat steeper Furthermore, in Sect. 5 we emphasised that 
the uppermost values in Fig.|6]practically represent a ceiling, as 
a consequence of the almost flat progression of the C with in- 
cresing 2 beyond about 300 MqPc^^, irrespective of Mbh (see 
Fig- El- Thus, in the context of our model, there are limits on the 
values that fabs can achieve, no matter the assumptions we are 
willing to make. 

As an exercise, let us assume that, in the local universe, the 
S distribution adopted in the previous two sections applies only 
to the late type galaxies, the early types being comparatively, on 
average, much poorer in molecular gas (Sect. 4). For simplicity, 
we assume that E is practically zero in the early type galaxies. 
Because the SMBH above 10^^"** Mq are mainly hosted by this 
type of galaxies, above Lx=10"^-''^ this assumption implies a ma- 
jor drop in fjbs, as illustrated in Fig. [12] at z=0.1, and the dis- 
crepancy with respect to the best fit by LF2005 is much stronger 
than in Fig. |9] Thence, to obtain a prediction closer in the local 
universe, as well as with increasing z, to the results in LF2005, 
one would need to assume that the presence of obscuring gas in 
at least the more massive early type galaxies (with the appropri- 
ate values of E and organised as the model requires) is correlated 
with the nuclear activity and E itself increases with z- One might 
adopt as a hint, in favour of the latter assumption, the increase 
with z in the "volume-averaged" star formation rate (SFR) (see 
also H2004 for further considerations). However, this is not so 
straightforward, because the evolution of the X-ray luminosity 
function of AGN is not only in their luminosity but also in their 
number density: thus, what occurs in the unit volume does not 
necessarily apply to the individual objects. Perhaps the simplest 
speculation, that one can naively propose, is that the quantity E 
effectively present, or at least surviving during the activity phase. 
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increases with A. However, as noted above, it would remain hard 
for our model to explain values of fats greater than about 0.3 at 
the highest luminosities. 

Unfortunately, no theoretical picture of the gas spatial distri- 
bution and kinematics, the interplay with the SMBH growth, on 
the one hand, and the SFR on the other, has been developed so 
far to a degree of (physically consistent) detail necessary to sen- 
sitively adjust the prediction of fabs- On the observational side, 
we lack sufficient data to proceed empirically, as we have done 
with A. 

We finally need to point out the non-negligible difference 
between the optical and the X-ray selection of AGN. Following 
Rigby et al. (2006) (see also Fiore et al. 2003, Cocchia et al. 
2007 and references in all three papers), we stress that a sub- 
stantial fraction of AGN selected through their X-ray emission 
and classified as such on the basis of their Lx would be missed 
through an optical spectroscopic classification (this fraction in- 
creases from about 15% locally to about 50% at z around 1). 
These objects could perhaps reduce the difference between our 
predictions and the fahsiz) found by LF2005. Such a possibility 
could be tested only after a properly complete investigation of 
the A of these objects is carried out. 
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Fig. 12. Same as Fig. |9] with the addition of the fabs (full line 
histogram) predicted assuming 2=0 in early type galaxies. 



8. Conclusions 

In the light of recent new results on the behaviour of the 
Eddington ratio /I as a function of the SMBH mass and of the 
redshift, we have revisited the predictions of the model pre- 
sented in Papl. This model naturally predicts an anti-correlation 
between the fraction fabs of X-ray obscured, Compton-thin 
AGN, and Mbh, provided there is interstellar molecular gas with 
appropriate values of surface density 2 approximately within 
25 to 450 pc (depending on Mbh) from the active nucleus. 
We used the newly found A{Mbh,z) dependence (z from the 
local universe to 0.7) and a distribution of 2, irrespective of 
the morphological type of the host galaxies and invariant with 
z, to convert that prediction into one versus Lx- The results 
appear to reproduce, at least qualitatively, the anti-correlation 
between fabs and Lx claimed by several authors, but not the 
still controversial claim of an increase in fabs with z- The latter 
increase, particularly in the form claimed by LF2005 on which 



we did concentrate, might possibly be reproduced by assuming 
that the typical values of E do increase with the redshift and, 
furthermore, that in the early type galaxies, at least in the more 
massive ones, the presence of this quantity is correlated in time 
with the activity of thek nucleus. At the highest luminosities, 
however, according to our model, fabs can hardly exceed 0.3. 
For the sake of completeness, it should be noted that the results 
just summarised are sensitive to the presence of a cut-off in the 
A(Mbh,z) distribution, and they correspond to the case where 
the cut-off is placed at the physical limit A=l. 

Steps forward, that are relevant not only to the model dis- 
cussed in this paper can be made at least in two ways. The first 
and most obvious one requires systematic investigation of the 
amount of molecular gas in an appropriately selected sample of 
early type galaxies, both active and inactive, out to redhifts of 
0.5 at least. The line strength and profile may not tell us enough 
about the spatial distribution, but much about the amount and the 
kinematics of the gas. The former property needs to wait for ob- 
servations with a large interferometric array of telescope, such 
as ALMA. The second step is somewhat indirect and therefore 
requires interpretation. The method is described and applied to a 
handful of broad-line objects in Maiolino et al. (2007). It allows 
an estimate of C to be inferred for each object on the basis of 
measurements of infrared spectral features due to dust reradia- 
tion of the ultraviolet from the active nucleus. It is worth noting 
that Maiolino et al. (2007) find an anti-correlation between C and 
either the luminosity or the mass, thus leaving completely open 
at present which of the two parameters is the one "physically" 
driving the anti-correlation. This sort of degeneracy could be re- 
solved by comparing the covering factor of objects with similar 
Mbh and widely different values of Lb, hence of A. According to 
our model a similar value of C should be found in such objects. 
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